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1 Introduction 

Let A he d X n matrix with entries in Nq = {0,1,2,...}. The i-th column of 
the matrix A is denoted by m. For a given /3 £ Nori = Vitti \ Vi £ No}, 

we consider the set of points 

j:-(/3) = {u G I riii =/?}. 

We are interested in the probability distribution /Z{p-,p) on F where u £ T, 
p G R" and Z is the normalizing constant defined by 

n 

uGJ- k—1 

When Qi lies on a hyperplane h{y) = 1, the probability distribution can be 
regarded as a conditional probability distribution induced from the multinomial 
distribution Y^'k=iPi — 1- We will call this probability distribution A- 

distribution following the celebrated work by P.Diaconis and B.Sturmfels on 
constructions of Markov bases of F in terms of Grobner bases of the affine toric 
ideal defined by the matrix A m. 0. The class of h-distributions includes 
hypergeometric distributions on contingency tables with a fixed marginal sums. 

We will discuss on an exact evaluation of the normalizing constant Z and 
its derivatives in this paper. The normalizing constant Z satishes the A- 
hypergeometric system, which is substantially introduced by Gel’fand, Kapra¬ 
nov, Zelevinsky in the late 1980’s (see, e.g., m), and contiguity relations (recur¬ 
rence relations) such as ^{j5\p) = Z{l3 — ai,p). By the system and contiguity 
relations, an discrete analogy of the holonomic gradient method (HGM) intro¬ 
duced in m will give an efficient and exact method to evaluate the normalizing 
constant and its derivatives. For the matrix A standing for the 2 x n contingency 
tables, the maximal likelihood estimate (MLE) by the HGM is discussed in the 
recent thesis by M.Ogawa m- The MLE needs evaluations of the normaliz¬ 
ing constant and its derivatives, which is done by contiguity relations for 2 x n 
contingency tables given in [3] . 
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In this paper, we will discuss on the discrete HGM for a general matrix 
A. The key step for it is a construction of the Pfaffian system for the A- 
hypergeometric system. We utilize the result of [5] , Macaulay type matrices m. 
and the Hilbert driven Buchberger algorithm HZ] for an efficient construction. 
Timing data are given to compare several methods. 


2 Macaulay Type Matrices 


The discussion of this section is well-known, but we do not find a relevant 
literature for our application. Then, we give a brief explanation on Macaulay 
type matrices (which has an origin in [ini P-8]). Let P = K[xi,... ,a;„] be the 
ring of polynomials of n variables where iti is a coefficient field. We suppose 
that a zero-dimensional ideal / of P is generated by polynomials /i,..., /m- Let 
g be an element of I. Then, g can be written as 

g = Cydy 1 fj = Cijtijfj (1) 

j \ i / ij 


where Cy G K and tij is a monomial. Let S' be a basis, consisting monomials, 
of P/I as a vector space over K. S is a finite set. We assume that a base S of 
P/I is given. In fact, in our application to H-hypergeometric systems, the base 
S is determined by an algorithmic method in [5] without computing a Grobner 
basis for a given H-hypergeometric ideal. We define M the union with respect 
to i's and j’s of the set of the monomials which appear in tijfj. Set M = M\S. 
Let us construct a matrix F of which columns are indexed by M U S' and whose 
each row stands for the coefficients of tijfj. We denote the row consisting of 
the coefficients of tijfj by C{tijfj). For g in dTj), C{g) is the row vector of 
the coefficients of the polynomial g, which is simplified, indexed by M U S. It 
follows from the relation m that the vector C{g) can be written as a linear 
combination of the rows of the matrix F. 

Let us present an example to illustrate our definitions. Put /i = X 1 +X 2 — 4, 
/2 = X 1 X 2 — 1. Let < the lexicograhic order such that xi > X 2 . Then, the set of 
the standard monomials S is {xf, ^ 2 , X 2 ,1}. Put g = X 2/1 —xi /2 = Xi-l-xf — 4 x 2 . 
Then, M = {xfx 2 ,Xi}. we have 


and 


MUS 

xfx2 

Xi 

xl 

xl X2 

C{x 2 fl) ( 

1 

0 

1 

0 -4 

C{xif2) 

1 

-I 

0 

0 0 

C{g) V 

0 

1 

1 

0 -4 

MU5 

x\x2 

Xi 

xl 

x| X2 

C{x 2 fl) 

( 1 

0 

1 

0 -4 

C(xi/2) 

V 1 

-I 

0 

0 0 


We suppose that the support of g, which is an element of I expressed as 
o, is S and a monomial t in M. We suppose that the coefficient of t in 5 
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is normalized to be 1. The monomial t can be uniquely written as a linear 
combination of the elements of S over K modulo / as t = mod/, 

Cs € K. Then we have g = t — X^sgS ^ is in the row span of 

F, we can obtain C{g) by constructing reduced row echelon form of the matrix 
F of C{tijfjys with the index order M > S. In fact, let e be the number 
of echelons (the number of the leading terms of the reduced Grobner basis of 
the system of linear equations defined by the matrix F). Then, the rank of F 
is equal to e. If C{g) is not contained in the reduced row echelon form, then 


the rank of 


is e + 1. It contradicts to the assumption that C{g) is 


F 
C{g) 

contained in the row span of F. 

For polynomials a positive number T, and a basis S of P/I 

consisting of monomials, we construct a matrix by row vectors of the form 
i = 1,... ,m where t runs over the set of the monomials of which degrees 
are less than or equal to T. We will call the matrix the Macaulay type matrix 
of degree T for /i,..., /^ and S. We denote the matrix by Ft- 


3 Macaulay Type Matrices to Derive Pfaffian 
Systems 

The method finding C{g) by computing the reduced row echelon form from the 
Macaulay type matrix has no advantage to Grobner basis methods without new 
ideas as in the Fa algorithm and its related algorithms, however, as we will 
see in this paper, this method has an advantage for the numerical analysis of 
^-hypergeometric systems in the ring of differential operators or in the ring of 
difference-differential operators without them. Consider the ring of differential 
operators with rational function coefficients R' = K'{di,... ,dn) where K' = 
C(xi,..., Xn), which is abbreviated as C(a;). Let / be a zero-dimensinal left 
ideal in R', S a basis of R/I over K'. Let s be an element of S. Then, 
diS is uniquely written as the linear combination of S over K' modulo I as 
diS = X^tGS mod /, Ct € K. Let us denote by g the expression diS — 

We call g the Pfaffian operator for di, s and S. We want to find all Pfaffian 
operators for all pairs of di, i = 1,... ,n and s G S. We utilize Macaulay type 
matrices explained in the section [21 which can be defined analogously in R'. 
Let F = Ft be the Macaulay type matrix constructed from t^-’s whose total 
degrees are less than or equal to T and fj's which are generators of I. When the 
total degree T is sufficiently large, then the ideal element g, of which support 
is in S' U {Sis}, can be obtained by computing the reduced row echelon form 
of F with an index order such that M > S and diS is the smallest index in 
M. We note that the linear algebra elimination over K' to obtain the reduced 
row echelon form does not perform differentiations of elements of K' = C(a;) 
and then the numerical coefficients of g when the variable x is specialized to 
a generic constant vector X can be obtained by specializing the matrix F to 
X = X and by the linear algebra elimination. This will be the key idea to apply 
Macaulay type matrices to evaluate A-hypergeometric polynomials, in which we 
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avoid a Grobner basis computation. 

Let A he & d x n matrix. The i-th column of A is denoted by at and the 
{i,j) element of A is denoted by aij. We assume that G Ng and the rank of 
A is d. Let us consider Macaulay type matrices when I is the ^-hypergeometric 
ideal Ha{c). We regard the parameters c^’s as indeterminates, and we put 

K = C{ci,...,Cd,xi,...,Xn) and R = K{di,... ,dn). 

Generators of Ha (c) consist of the Euler operators Ej — Cj = J2k=i ^jkXkdk — Cj , 
j = l,...,d and the toric ideal I a = (9" — d'"\Au = Av,u,v G Ng). It 
is important to distiguish the left ideal RHa{c) and the left ideal R'Ha{/3), 
jd G The former ideal lies in R where Cj’s are indeterminates and the 

latter ideal lies in R'. We consider a Macaulay type matrix for Ha{c). We will 
construct a smaller matrix F' than the matrix F to obtain Pfaffian operators 
in the following way. 

For a given term order <, the reduced Grobner basis of I a with respect 
to the order < is denoted by G. A basis of R/RHa{c) can be found by the 
algorithm in [5] from G. We denote it by S. Take a monomial s from the set 
S'. 

Algorithm 1 

Input: a matrix A which defines an A-hypergeometric system and a sufficiently 
large integer T. 

Output: A matrix F^ with entries in K , which will be called the Macaulay type 
matrix of degree T of the A-hypergeometric system. 

1. Gompute a Grobner basis G of Ia- 

2. Find a basis S by the algorithm in [5]. (A probabilistic version of this 
method in [7] can be used and is more efficient.) 

3. Let Nt be the set of the monomials in dfs whose total order is less than 
or equal to T. 

4. Multiply an element t of Nt to an Euler operators Ej — Cj in the ring of 
differential operators R. Reduce the element t{Ej — Cj) by G and obtain 
the remainder. By the remainders from all the pairs from Nt and the 
Euler operators, we construct the matrix F' as in the section [7] and set it 
F^. 

Theorem 1 1. When T is sufficiently large, the reduced row echelon form of 

F' = Ff contains the Pfaffian operator diS — irreducible 

by G. 

2. Let f be a solution of Ha{c). The numerical value {dis) • f at a generic 
point X = X G Q", c = /3 G can be obtained from the numerical 
values of s • f, s € S at a point x = X,c = (3 by computing the reduced 
row echelon form of the numerical matrix F'|a;=x,c=/ 3 - 
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We note that when diS is reducible by G, it is equal to an irreducible monomial 
9" modulo G. When T is sufficiently large d"" — reduced 

row echelon form of F'. Then, we can obtain the Pfaffian operator diS—'^^^g Ctt 
from the reduced row echelon form. 

Proof. 1. The set S is a basis of R/RHa{c) over K = C{ci,... ,Cd,xi,... ,Xn). 
We multiply Nt from the left to the generatros Ej — Cj and 9“ — O'" of Ha{c). 
M is the set of the momonials in the obtained elements of R from Nt ■ {Ej — Cj) 
and Nt • G. Put M = M\S as va the sectionj^l The set G is a Grobner basis of 
the toric ideal I a- We devide M into two disjoint sets Mr and its complement 
Mf where the set Mr is the set of reducible monomials by G. We define the set 
Mi the union of Mf and the set of irreducible monomials obtained by reducing 
elements of Mr by G, which is denoted by Red(Mr), minus the set S. Thus, for 
9“ G Mr, there exists 9" G Mi U S such that 

9“ - 9" G I A. (2) 

The subset of Mi , of which element appears in the column index of E ', is denoted 
by M'. The set M U (Red(Mr) \ S) consists of the disjoint sets Mr, Mi \ M', 
M'. Construct the Macaulay type matrix E from these binomials 0 and the 
obtained elements by the multiplication of Nt to the generators of Ha{c). We 
sort the columns of the matrix by the order Mr > {Mi\M') > M' > S. When T 
is sufficiently large, the reduced row echelon form of the Macaulay type matrix 
E contains the Pfaffian operator diS — Gt as we have seen in the previous 

section. The reduced row echelon form of the smaller matrix F' is contained in 
the reduced row echelon form of F. Since diS G Mi n M' for sufficiently large 
T, the Pfaffian operator is contained in the reduced row echelon form of F'. 

2. The steps to find the reduced row echelon form do not have steps of 
multiplication in the ring of differential operators R. Then, a specialization of 
X and c to numbers and constructing the reduced row echelon form commute. 
Q.E.D. 

Example 1 Consider the A-hypergeometric ideal generated by 


a;i9i + X2d2 + x^d^ + X/^d/i — ci, X2d2 + 2^494 — C2, x^d^ + 2:494 — C3, 
828^ - 8 i 8 a . 

For the reverse lexicographic order such that 9 i >82 >83 >84, the basis 
by the algorithm is S' = { 1 ,94}. Put the degree T = 1 . We multiply the 
monomials in Nt = { 1 , 9 i, 92 , 93 ,94} to the generators. The table below is the 
result of the multiplication where the index iii2 ■ ■ ■ im in the top row stands for 
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the monomial index 0 denotes the monomial 1. 


M' Mr M' S 

11 12 13 14 22 23 24 33 34 44 1 2 3 4 0 



We put 


Mt = {23,114,123,124,134,144,223,233,234} 

M' = {1,2,3,11,12,13,14,22,24,33,34,44}. 


Note that Mt = U (M* \ M') in the proof. The join of the two tables is 
the matrix F and M = Mt U M'. Let us use the order of indices such that 
Mt > M' > S. Apply the Gaussian elimination with this order to the first 
table. We eliminate elements of the column standing for the index 23 by the 
last row and then remove the last row. Then, we obtain the following matrix 
which agrees with the matrix F^, (T =1) in the algorithm. 



We can see, by a calculation, that this matrix can be transformed into the 
reduced row echelon form whose rank is 12. The reduced row echelon form 
contains the Pfafhan operators. 


4 Evaluation of 74-Hypergeometric Polynomials 
by Macaulay type matrices 

Let R' = K'{di ,..., dn), K' = Q{xi ,..., Xn) be the ring of differential operators 
with rational function coefficients. Let I be a zero dimensional left ideal of R' 
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and S = {si,...,Sr} a basis of R/I as a vector space over K'. We assume 
that the set S consists of monomials of d. If no confusion arises, the column 
vector (si,..., is also denoted by S. There exists a matrix Pi with entries 
in K' such that diS — PiS = 0 modulo I. The system of differential equations 
diY — PiY = 0, i = 1,... ,r where T is a column vector of unknown functions 
is called a Pfaffian system. A basis S and the matrix Pi can be obtained by 
computing a Grobner basis of I in the ring of differential operators R' (see, e.g., 
[H Chapter 6]). Let Ha{/3) be the A-hypergeometric ideal. Let Z{l3-,x) be an 
A-hypergeometric series for the matrix A and a generic parameter vector (3. It 
is well known that a contiguity relation 

di • Z{l3;x) = Z{f3 - af,x) 

holds under a suitable normalization of the hypergeometric series, see, e.g., 
|151 p.xy]. In particular, the relation holds when Z is the A-hypergeometric 
polynomial 

1 n n 

E u\ = X{u.\,x^=X{xT (3) 

Alii—/3 ,uGNq i—1 i—1 

for /3 £ No A = Nofli- It follows from the contiguity relation that we also 
have a contiguity relation for derivatives of Z 

di • {sj • Z){P\x) = {sj • Z){I3 - ai]x) 

We define the column vector F by (s • Z | s S A). When 5” is a basis for both 
R'/R'Ha{I3) and R'/R'Ha{I 3 — cii), it follows from the contiguity relation and 
the Pfaffian system, we have the following identity 

F(/3 - ai;x) = di •Y{j3;x) = Pi{P]x)Y{l3;x), (4) 

which is called a recurrence relation or a difference Pfaffian system for the vector 
valued function Y to the direction a^. 

We suppose that the value F(c, a;) at c = /3 S NqA and x = X £ Q" is given. 
Once we construct the difference Pfaffian system, the value of F at c = /3 + 
and X = X is obtained by 


P,{l3 + a„x)-^Y{f3,X) (5) 

when the inverse matrix of Pi exists. When c = 0, the hypergeometric poly¬ 
nomial Z{0;x) is equal to 1. Then, the vector valued function F(0;x) is 
(I, 0,..., 0)^. We apply the recurrence relation ([5]) from F(0;x) iteratively 
for i’s. Then we can obtain the exact value of Y{J2m=i ^ Nq as long 

as the inverse matrices of Pi’s exist and S' is a basis for all R'/R'PIA{d). This 
method to evaluate hypergeometric polynomials is called (difference) holonomic 
gradient method ((difference) HGM) as an analogy of evaluating normalizing 
constants for unnormalized probability distributions by utilizing holonomic sys¬ 
tems of differential equations m. ID- 
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Finding a recurrence relation of Z for indeterminates xi,..., and ci,..., Cd 
is possible by deriving a PfafHan system by a Grobner basis, but it requires huge 
computer resouces. When we want to evaluate Y (c; x) at a vector of rational 
numbers x = X on a line c = l3 + kH parametrized with a parameter k and a di¬ 
rection vector H from the value Y{j3]X), we only need the matrix Pi restricted 
on X = X and c = P + kH. In order find such restricted Pi, we can apply 
the method of Macaulay type matrix and avoid Grobner basis computation for 
indeterminates x and c. This idea is illustrated in the Figured] 


Restriction of P on a line c = P + kH and x = X 


Method of Maibaulay type matrices 


Grobner basis method 


Specialization of c to /3 -f kH, x io X 


ideal I 


PfaSian system P oi I 


Figure 1: Restriction to c = fi + kH and x = X 


Algorithm 2 (Method of Macaulay type matrices for A-hypergeometric sys¬ 
tems ) 

Input: a matrix A, a starting point of the parameter vector /3 £ Nq, a vector 
of numbers X G Q”, a basis S of Rn/RnHA{c) consisting of monomials of d, a 
direction H G Nq such that H G ^oOi- 

Output: the recurrence relation 

Y(k - 1) = R(k)Y(k) (6) 

with respect to k for the vector valued hypergeometric functions Y (fc) = (S' • 
Z){l3 + kH-,X). 

1. Find hi G Nq such that H = (express the direction H in terms 

of Oi). This step can be performed by solving the integer program problem 
of minimizing w ■ h under the constraint Ah = H, h G Nq for a weight 
vector w. 

2 . Let Gij^ be the Grobner basis of the toric ideal I a in the Algorithm dl Let 
Sj be the normal form N{d^Sj,Gi^) of d^Sj with respect to the Grobner 
basis . Ghoose sufficiently large number T such that all s', j = I,..., r 
are contained in M' U S in the Algorithm [T| 

3. Gonstruct the Macaulay type matrix by the Algorithm [T| 







4. Restrict F!p to x = X and c = ^ + kH and set the obtained matrix If 
the reduced row echelon form of contains the element standing for 

r 

= 'y ] ^3 1 ^ij S Q(fc) (7) 

i=i 

for all j = 1,... ,r, then go to the step 5 else increase T and go to the 
step 3. 

5. Define the matrix R{k) by (ty). 

/ 1 1 1 1 \ 

Example 2 Put A = I 0 1 0 1 and set the direction H = (1,1,1)^, 

\ 0 0 1 1 / 

which stands for 94 , put the starting point f3 = (3,2,1)”^, and put X = 

(1,1,1/2,1). As we have seen in Example[I] we can choose S = (1, . Then, 

we obtain 

( Z{P-{k-l)H-,X) \ 

\ {d^.Z)[P-{k-l)H-X) ) 

^ ( 0 1 V Z{f3-kH-,X) \ 

\ -2{k + l){k + 2) 3k + 5 J \ {d4» Z){P - kH;X) J 

by our algorithm. This recurrence relation holds for any hypergeometric poly¬ 
nomial for the parameter c = P + kH. The hypergeometric polynomial Z{c',x) 

k + 1 

is -f 2^2*32^4 for C = (3 -b fc, 2 -f fc, 1 -f k)'^. 

5 Hilbert Driven Algorithm for A-Hypergeometric 
Ideals 

In the previous section, we discussed a method to evaluate hypergeometric poly¬ 
nomial by contiguity relations and Macaulay type matrices. The evaluation by 
contiguity relations (the difference HGM) can also be performed if we have a 
Grobner basis for the ideal Ha{c). The rank of RHa{c) where Cj’s are inde- 
terminates is equal to the normalized volume of A by Adolphson’s theorem [ 1 ] , 
then the Hilbert driven algorithm by Traverso m can be applied to the A- 
hypergeometric ideal Ha{c). Let us summarize this method as an algorithm. 


Algorithm 3 (Hilbert driven method for A-hypergeometric systems) 

Input: a matrix A standing for the A-hypergeometric system and (3 € NqA. 
Output: recurrence relations (difference Pfafhan systems). 

1. Evaluate the normalized volume of A. Put it r. 0 

^The number r can be obtained by evaluating the multiplicity of I a or by evaluating the 
volume by geometric methods. 
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2. Compute the Grobner basis G of Ha{c) in the ring of differential opera¬ 
tors R — Q(ci ,... ,Cd,xi,..., Xn){di ,..., dn)- We use the Hilbert driven 
algorithm to avoid unnecessary 5'-pair checks. In other words, we stop 
the Buchberger algorithm once the number of standard monomials of an 
intermediate Grobner basis equals to r. 

3. Find a path from /3 to /3' which is near to 0 by the Algorithm 0] 

4. Let S be the vector of the standard monomials for G. Compute matrices 
Pi such that diS — Pi S modulo G by the normal form algorithm only for 
the indexes i which appear in the path. 

5. We have recurrences Pi{c + ai)Y (c + ai) = Y(c) on the path. 

At each step of applying recurrences, we specialize parameters c and x in Pi to 
numbers. It is our heuristic observation that Ci should not be speicialized to 
numbers or parametric polynomials during the Grobner basis computation and 
the normal form computation. It seems to make the computation slower. 


6 Comparison of Methods to Evalute Hyperge¬ 
ometric Polynomials 


In the latter two sections, we have illustrated two methods to evaluate A- 
hypergeometric polynomials numerically. We will compare the two methods 
and the method by enumerating the fiber {it € Nq | Au = /3}. 

The advantage of the method of Macaulay type matrices is that the step 4 
of the computation of the reduced row echelon form is performed in the held 
Q(fc). On the other hand, the Grobner basis method needs computation in the 
ring Q(c, x)(di ,..., dn), which contains d + 2n indeterminates. 

All timing data in this section is taken on a machine with Intel Xeon CPU 
(2.70GHz) with 256 G memory. 


Example 3 

Risa/Asir. 


Timing data of test_clllc_conti2() in the package 


for 


Put A = 


/3 = (3,2,1,1)^ 


/ 1 1 1 1 1 1 1 1 \ 

0 10 0 110 1 

0 0 10 10 11 

V 0 0 0 1 0 1 1 1 y 

A = (1,1/2,1/3,2/3,1,1,1,1)5= (1, as, de, dr, dg, dif. 

Our benchmark problem is to get the value of S • Z at x = X and c = 
d + k(3, 1,1,1)^. The contiguity for (1,1,1,1)^ is applied to obtain the value 
at (3 -\- fc(l, 1,1,1)^ and next the contiguity for (1,0, 0,0)^ is applied to obtain 
the value at j5 -\- fc(3,1,1,1)^ from the value at (3 -\- fc(l, 1,1,1)^. 
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THe following table is timing data of the method of Macaulay type matrices. 


k 

Time (second) 

0 

1.45 

10 

1.48 

20 

1.70 

30 

1.79 

40 

1.89 

50 

1.99 

60 

2.16 

70 

2.35 

80 

2.58 

90 

2.82 

100 

3.17 


The timing for the case k = 0 is the time to construct recurrence relations for 
the direction (1,1,1,!)^ and (1,0,0, 0)^. When k = 10, our implementation 
outputs 


Val=[30318066527332447242457/89619251224349337722522492794306560000, 

. . .]. 

From the output, we obtain, e.g., values 

Z = 30318066527332447242457/89619251224349337722522492794306560000 
= 3.38... X 10-1® 

Eps] = ds* ZjZ = 52047189429143224956864/30318066527332447242457 

= 1.71... 

where E[Us] is the expectation of the random variable Us G N® which satisfies 
AU = /3 + k{3, 1,1,1)^. We compare these data with other two methods for 
exact evaluation. First one is an exhaustive enumeration of 17 G Nq satisfying 
AU = (3 + fc(3,1,1,1)^. It can be easily enumerated for this by a nested loop 
of four for statements, because when non-negative integers U5,Uq,Uy,Us are 
given, other 67 ’s are determined uniquely as integers and the admissible range 
of Us,... ,Us can be described in terms of k. Moreover coefhcients can easily be 
evaluated by recurrences. 
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Figure 2: black square=Macaulay type matrix, white square=Hilbert driven 
GB, white circle=series, horizontal=A:/10, vertical=seconds 


k 

Time (second) 

number of fibers 

0 

0 

5 

10 

0.0200 

1946 

20 

0.132 

18436 

30 

0.476 

76976 

40 

1.42 

220066 

50 

3.70 

505206 

60 

7.68 

1004896 

70 

15.9 

1806636 

80 

30.9 

3012926 

90 

56.0 

4741266 

100 

88.7 

7124156 


The last table is timing data of evaluating the normalizing constant in the 
benchmark problem by recurrence relations derived from Pfaliian operators ob¬ 
tained by the algorithm |3l The Hilbert driven Grobner basis computation in 
the step 2 takes 5.50 seconds on the yang package for the Risa/Asir. If we do 
not use the Hilbert driven method in the step 2, it takes 2069.78 seconds on the 
yang package. In our benchmark problem, we may only obtain Pi and Pg- It 
takes 6.80 seconds to obtain them in the step 3. 
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k 

Time Recc (second) 

Time Recc + GB (second) 

0 

0.02 

12.31 

10 

0.06 

12.35 

20 

0.072 

12.37 

30 

0.14 

12.43 

40 

0.27 

12.56 

50 

0.34 

12.64 

60 

0.49 

12.78 

70 

0.62 

12.92 

80 

0.83 

13.13 

90 

1.05 

13.34 

100 

1.32 

13.62 


In the table, “Time Recc” is the time to apply the recurrence relations to obtain 
the value of the normalizing constant. “Time Recc + GB” is the sum of the 
“Time Recc” and the time of the algorithm |3l 

The Figure [2] tells us a strategy to evaluate A-hypergeometric polynomials 
numerically. When j3 is small, then the enumeration of the fibers will be the 
best method. When [3 is large, then the method of Macaulay type matrix will be 
the best choice. When we need to evaluate the hypergeometric polynomial for 
several /3’s, then the method of Grobner basis by the Hilbert driven algorithm 
will be the best choice. 


7 Finding a Path to Apply for Recurrence Re¬ 
lations 

Suppose that /3 G NgA. We want to find a path in NqA from /3 to jS' closer to 

0 

Algorithm 4 

Input: A and /3 € NqA. 

Output: a path to /3' G Nq 

1. Express /3 as /3 = niat^ + ■ • • + Ukai^, rii G Z>o by reducing by a 
Grobner basis of ■ 

2. Find the maximum in the set 

{m I path /3 - 7u -t /3 - 7u - mai. 

lies in NqA for all 7« = X] ^kak, G S and m < nj} 

When the set is empty or ,9 = 0, return the path and /? as the terminating 
point P'. 

3. Add the pair to the path. Put rij = rij — m and /3 = /3 — maij . Go 

to the step 2. 
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The definition that “path /3— 7 „ ^ /3— 7 „—mai^. lies in NqA” is that S 

No^, /3 - 7 „ - e NqA, ..., /3 - 7 ^ - ma^- e NqA. 

Let S' be a basis of R/RHa{c) consisting of monomials in di's and diS — PiS, 
i = 1,... ,n is the Pfaffian operator for S. The singularity polynomial of the 
Pfaffian operator is the least common multiple of the denominator polynomials 
of the elements of the Pi’s. It is a polynomial in C[c, a;]. The basis S is called 
a good basis for NqA when the following two conditions are satisfied. 

1. The singularity polynomial for the Pfaffian operator does not contain the 
variables Ci’s. 

2 . S is a basis of R'/R'Ha{I 3), R' = C(x)(di,..., 9„) for all /3 G Nq^. 

Note that when A admits a rank jumping parameter /3 G NqA [3], there exists 
no good basis. 

The following theorem shows that the algorithm |4] and HGM reduce the 
evaluation of Z(/3; x) to Z{I3'] x) where /3' closer to 0 when the condition of the 
theorem holds. 

Theorem 2 When A is normal and S is a good basis, then the matrix Pi (/3 + 
af, x) has the inverse for j3 satisfying /3 — X) € NqA for all 9“ G S' and for 
any x out of a measure zero set. 

Proof. We will construct the inverse matrix of Pi. We denote by Ma{P) 
the left P„-module Dn/DnHAiP). Since A is normal, it follows from Saito’s 
isomorphism m that Ma{I3) and Ma{P + Oi) are isomorphic as the left £)„- 
module when /3 G Ng^. The isomorphism is given by 

Ma{I3) 9 f I—>■ £di G Ma{I3 + a^). 

In general, they are isomorphic when ft + Oi ^ V{Bi), which is the zero set of 
the b-ideal for the direction i [14]. When A is normal, we have (3 + Oi ^ V{Bi) 
for all i and /3 G NgA. Note that we have the induced isomorphism 

Hom_D^(MA(/3 + ai),0) 3 f 1 -^ di • f e Hom_D„ (M^(/3), O). 

Let 9“ be an element of S. When 9 — X) '^kO-i G NgA, there exists an operator 
£iu G Dn such that the following diagram of left P^-modules commutes and all 
arrows are isomorphisms. 

^a(/3 - X Ukak) < Ma(/ 3 + Ui - X Ukak) 
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Since the diagram commutes, we have 


= 0 inM^(/3 + ai), 


which implies iiud'^di — 9“ G DnHA{P + Ui). Let / be a local solution of 
Ma{13), which can be regarded as an element of Homu^ (M^(/3), O). Con¬ 
sider the column vector l/(/3) = (s • /1 s G 5'). Take a local solution g in 
Hom£)„(M^(/3 -I- ai),0). We have iiud^di • g — • g = 0. The function 

f = di • g is a. solution of Ma{I3)- Then, we have Yg{P + m) = (9“ *519“ G 
S) = • /1 9” G S). Suppose that = '^Cy{x)d'". Since S' is a 

good basis, d"" can be expressed as a linear combination of S with coefficients 
in C(x) modulo R'Ha{/3). Therefore, there exists a matrix Qi with ratio¬ 
nal function entries satisfying Yg^P + Ui) = Qi{x)Yf{j5). We note that Qi 
does not depend on the choice of the function g. It follows from the conti¬ 
guity relation that we have Yf{j3) = Pi{/3 -I- ai,x)Yg{P -|- Oi). Hence, we have 
Yg(/3 + Qi) = Qi{x)Pi{j3 + ai,x)Yg{P + ai). Make g run over a solution basis. 
Then, we can see that Qj{x)Pi{/3 + ai,x) = E out of an analytic set in the 
x-space. The analytic set defined by a non-zero holomorphic function is a mea¬ 
sure zero set. A countable union of measure zero sets is a measure zero set. 
Q.E.D. 

Note that we do not have an example of S which is not good when A is 
normal. It is an interesting question to find such example or to prove the 
existance of a good basis. 

8 Reduction to Non-Negative a^’s 

We have supposed G Nq in the previous sections. We assume that at G Z‘^ 
and there exists a linear form h G R" such that h{ai) = 1 for all *’s in this 
section. We will explain a method to reduce this case to the case ai G Nq. 

We define the vector p G Ng by 


Pi = -minj(ay,0). 


When h{p) = —1, we choose a unit vector e which is not orthogonal to h and 
rewritep by p-|-e. Then, we can assume that h{p) ^ —1. Define the new matrix 
A! by (oi -I- p,..., a„ -I- p). The toric ideals for A and for A' agree, e.g., [3 
Lemma 1.5.10]. 

Let us discuss on a relation between solutions of the Ha and Ha’ ■ When / 
is a solution of Ha{I3), we have 


n 



Multiplying hi and adding with respect to i, we obtain 


n 


'^Xjdjf = h{l3)f 


15 


by utilizing h{aj) = 1. Adding Pk(J2'j=i — ^iP)Pkf to the fc-th equation 

Ekf = Pkf, obtain 

'^{akj +Pk)xjdjf = {Pk + h{P)pk)f 

Hence, the solntions of Ha{P) are solutions of Ha'{P') where /?' = /3 + h{P)c. 
In particular, the hypergeometric polynomial for Ha{P) is the hypergeometric 
polynomial for Ha'(P')- 
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